Well hello there, partner!
My name is Conor Higgins, and I’m a first year masters student in the Ecology program at UC Davis.
Generally speaking, my thesis work investigates community level interactions between ungulate herbivores and the associated plant community in an arctic ecosystem, and how these interactions may be impacted by climatic drivers. My field work is conducted at our lab’s long term study site near Kangerlussuaq, Greenland. I conduct near daily counts of the two large ungulate herbivores there (caribou; Rangifer tarandus, and muskoxen; Ovibos moschatus) as well as collect data on the timing of various plant phenological stages. My current research questions revolve around the reciprocal interactions between muskoxen and the plant community. More specifically:
and reciprocally
Muskoxen: image taken from https://news.nationalgeographic.com/content/dam/news/2018/01/18/muskoxen-climate/01-muskoxen-arctic-climate.ngsversion.1516271410305.adapt.1900.1.jpg
#Use the storms data (included in RStudio) to create a plotly graph of the relationship between wind and pressure, where the status of the storm is indicated by a color.
library(tidyverse)
library(plotly)
ggplotly(ggplot(data=storms)+
geom_point(aes(x = wind, y = pressure, color = status), pch=21, alpha= .05) +
scale_fill_viridis_c()+
theme_bw()+
labs(title = "Hurricane Wind Speed v Pressure", x = "wind speed"))
#Create a table that identifies the mean wind, pressure, ts_diameter, hu_diameter of each status of storm (remember to remove NAs!). Use the package htmlTable. Round each mean to only two decimal places (Hint look up the function round)
library(htmlTable)
storms <- storms
table <- storms %>%
drop_na(ts_diameter, hu_diameter) %>%
group_by(status) %>%
summarise(mean_wind = mean(wind), mean_pressure = mean(pressure), mean_ts = mean(ts_diameter), mean_hu = mean(hu_diameter))
?round
table$mean_wind <- round(table$mean_wind, 2)
table$mean_pressure <- round(table$mean_pressure, 2)
table$mean_ts <- round(table$mean_ts, 2)
table$mean_hu <- round(table$mean_hu, 2)
#CHALLENGE Find the duration, in number of days, of every hurricane from 2010 and later, and then use one of the map functions from purrr to write a sentence saying “Hurricane X lasted Y days” for each of these storms.
library(purrr)
stormduration <- storms %>%
filter(year >= 2010, status == "hurricane") %>%
group_by(name) %>%
summarise(duration = diff(range(day)))
map2_chr(stormduration$"name", stormduration$duration, function(x,y) paste("Hurricane", x, "lasted", y, "days"))
## [1] "Hurricane Alex lasted 29 days"
## [2] "Hurricane Arthur lasted 2 days"
## [3] "Hurricane Chris lasted 0 days"
## [4] "Hurricane Cristobal lasted 3 days"
## [5] "Hurricane Danielle lasted 7 days"
## [6] "Hurricane Danny lasted 2 days"
## [7] "Hurricane Edouard lasted 4 days"
## [8] "Hurricane Ernesto lasted 1 days"
## [9] "Hurricane Fay lasted 0 days"
## [10] "Hurricane Fred lasted 30 days"
## [11] "Hurricane Gonzalo lasted 6 days"
## [12] "Hurricane Gordon lasted 2 days"
## [13] "Hurricane Humberto lasted 2 days"
## [14] "Hurricane Igor lasted 9 days"
## [15] "Hurricane Ingrid lasted 2 days"
## [16] "Hurricane Isaac lasted 1 days"
## [17] "Hurricane Joaquin lasted 29 days"
## [18] "Hurricane Julia lasted 3 days"
## [19] "Hurricane Karl lasted 1 days"
## [20] "Hurricane Kate lasted 0 days"
## [21] "Hurricane Katia lasted 9 days"
## [22] "Hurricane Kirk lasted 30 days"
## [23] "Hurricane Leslie lasted 6 days"
## [24] "Hurricane Lisa lasted 1 days"
## [25] "Hurricane Maria lasted 1 days"
## [26] "Hurricane Michael lasted 5 days"
## [27] "Hurricane Nadine lasted 29 days"
## [28] "Hurricane Nate lasted 1 days"
## [29] "Hurricane Ophelia lasted 29 days"
## [30] "Hurricane Otto lasted 1 days"
## [31] "Hurricane Paula lasted 2 days"
## [32] "Hurricane Philippe lasted 4 days"
## [33] "Hurricane Rafael lasted 2 days"
## [34] "Hurricane Richard lasted 1 days"
## [35] "Hurricane Rina lasted 3 days"
## [36] "Hurricane Sandy lasted 5 days"
## [37] "Hurricane Shary lasted 0 days"
## [38] "Hurricane Tomas lasted 26 days"